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ON THE DESIGN OF STRUCTURAL MODELS TO STUDY 
THE THERMAL STRESS PHENOMENON* 


L. ALBERT SCIPIO II 


Department of Aeronautical Engineering 
University of Minnesota, Minneapolis, Minn. 


INTRODUCTION 


Elastic or structural similarity is generally very simple to formulate. 
However, the complexities associated with thermal stress, buckling and 
related behavior, plastic flow, use of material for the model different from 
that in the prototype, and problems associated with fluid-structural com- 
binations such as flutter, require very comprehensive analysis. 

A true structural model gives an accurate prediction of distortion as 
well as stress resulting from all four basic types of loadings: (1) axial, 
(2) shear, (3) torsion, and (4) bending. The structural characteristics of 
the prototype such as deflection, strains, stresses, buckling loads, and/or 
combinations thereof, could be obtained from a model that is geometri- 
cally similar, if all dimensions, connections, spacings, and so on have 
been reproduced to scale. 

In order to maintain elastic similarity, it is necessary that the mean 
deflection of the model and that of the prototype be similar and of the 
same scale N,.! In addition, it is required that the model simulate the 
buckling characteristic of the prototype and, in the case of wing struc- 
tures, the local buckling of the skin. 

Elastic similarity is discussed from the standpoint of two structural 
types: (1) solid and (2) wing-type structures. 

The following four basic assumptions are made: (1) the stress field is 
a homogeneous, isotropic, elastic medium, and creep is not present; 
(2) the elastic proportional limits of the material of the model or the 
prototype are not exceeded; (3) the stress fields are slowly varying func- 
tions of time; and (4) the physical properties of the material are not 
independent of temperature. 

The symbols used in this presentation are listed in TABLES 1 and 2. 


*The work reported in this paper was supported in part by Contract AF 
18-600-851 from the Air Research and Development Command, United States Air 


F orce,Wright-P atterson Air Force Base, Ohio. 
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TABLE 1 


ENGINEERING SYMBOLS 


A characteristic length, X,Y,Z Body forces 
thickness x,y,z Space coordinates 
Young’s modulus of elasticity a Coeificlent of thermal 
Component of total deformation expansion 
Component of external force re) Rib and web spacing 
Characteristic thickness € Component of deformation 
Scale factor, sectional forces y Component of shear deformation 
Outward normal to surface of kK Thermal diffusivity 
body oO Stress component 
Pressure i. T Shear stress component 
C°irections respectively —-«-51y_-—=s Components of unit invariant 
SUBSCRIPTS 
i,j Indices 
L Lower skin 
U Upper skin 
W Webs and ribs 


oe Ae = a | 


X,Y,Z Space coordinates 


TABLE 2 
AEROELASTIC SYMBOLS 


Young’s modulus 

Vertical translation, bending frequency effects 

Moment of inertia of wing segment about the elastic axis 
Mass of wing segment per foot 

Time 

Flutter speed 

Density 

Flutter frequency 


Subscripts 


Flutter effects a Torsional frequency effects 


1.0 SOLID STRUCTURES 
(1.1) Stress Field in a Solid Body 


Thermoelastic equations. If the temperature at the point p of the body 
is increased by 0, deformation can be produced without any corresponding 
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change of pressure. If forces are applied at the point in addition to the 
increase of temperature, the strain at the point consists of the thermal 
extension superimposed upon the strain connected with the stress. The 
thermal expansion, which all linear elements of the body undergo, is « 6, 
where « is the coefficient of thermal expansion of the body material and 
0 is the temperature increase. 

The stress-strain relation for a body strained by change of temperature, 
in addition to body forces, then, is given as follows: 


Dr ae) 
cee eve te +40 ye, 
1 2 (1 + v) 
sigh = lo, -—v(o,+0,)]+a0 Vers ry paar Tyz (ict) 
1 2(1 + rv) 
€s2 = [o, -—v (o, + oy)] +a Yzx = E Tzx 


General equations of equilibrium and compatibility. The general equa- 
tions of equilibrium of the surfacé stresses with the external forces on 
the surface can be written as: 


00x Oia are 
+ —_— —4+¥ 


= =F, 

ox dy dz 

Ody F: Ofzy rf Otyz = -F, (1.2) 
oy Ox dz 

Oo, OTzs Otyz = =r 


+ + 
Oz Ox oy 
where F,, F,, and F, are components of external force. 
The conditions of compatibility are: 


2 
a er 6 Dre ave at) pies 


Oy? . Ox? dxdy~— Ox Ox Oy Oz 


dey 07e, OaV ne C] ( OVy2 F, Wy 2 % ay 9 Oey 


Oz? ¥; dy? - Oydz ay ox oy dz 


07, 07¢, 07y,- O ( Vee - OVxz OVxy ‘ Oe, 


pp gz? dxdz dz ox oy dz Oxdy 
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Boundary conditions. The total on the surface of the body consists of 
normal aerodynamic pressure —P6ij and the viscosity stress of the gas 
(aerodynamic friction drag) p (dv;/dx; + Ov;/0x;). The boundary conditions 
are, therefore, 


a a 
be (— ut — \ Pa et (1.4) 
Ox; 


In addition, boundary conditions regarding the supporting structure of 
the prototype should be duplicated in the model. 


(1.2) Elastic Similarity Conditions 


On the basis of EquaTions 1.1, 1.2, and 1.3, we can write the follow- 
ing similarity equations for the stress field in a body produced by exter- 
nal forces and temperature: 


0; (a;;) ij e 0 x 
J =f [tw , ot ’ = , —. »* Oo, V; =] (1.5) 
E > Ee Dee L 


where e; = u, v, and w. Similarity equation 1.5 gives the following simi- 
larity conditions: 


Nx, = Ne; =N1 (1.6) 
Noi, = Neo, = N74 = NE =Np (1.7) 
VoNe=1 (1.8) 
Ny= 1 (1.9) 


2.0 WING STRUCTURES 
2.1 Mean Deflection of Wing Structures 


Consider an element of the wing structure (F1GURE 2.1) that is sub- 
jected to lateral forces and forces acting in the middle plane of the skin. 
Tangential forces are denoted by N,, N,, N,,, and N,x, and for the body 
_ forces or tangential forces per unit area of the median surface we have X 
and Y acting along the x and y directions respectively. 

By projecting these forces on the x and y axes we obtain the following 
equilibrium equations: 
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= —X 
Ox oy 
aN, ON, er, 
Ox oy 
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(2.1) 


The equation of equilibrium is formed by equating the summation of the 
normal and shearing force projections on the z axis to the load on the 


element. We can write that 


0*w 0*w 0?w ow [{ ON, 
p+, ta AN, ae; 2N xy : i 
Oy Oxdy Ox Ox 


ox? 


oe 
oy oy Ox 


+ 
Ox Ox? dy ody? 


ow { AN, a) _ OM aw ON, d?w , 


Substitution of EQUATION 2.1 gives 
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ON, 0?w ON, 0?w 
x + dy = 
Ox Ox? oy oy? 
(2.3) 
x o°w a 0?w Sy 0*w _x ow aye 
NP gat dy? cd oxdy ox oy 


EQUATION 2.3 represents the deflection of a plate, if there are forces 
in the middle plane in addition to lateral forces. The usual representation 


is 
En? ( ate 8 aie ir) 
12(1-v?) \ ox* i Ox*dy? as dy* a 
(2.4) 
2 2 a? Ow Ow 
ca O*w aii O*w + 2N,y w eee 
Ox? vs dy? Oxdy Ox Oy 


EQUATION 2.4 can be applied to wing structures which consist of two 
plates (an upper and a lower) connected by a common member or members 
with a spacing of 6 as shown in FiGuRE 2.2. 

We define V,, N,, Ny, =N,, and Ng as follows: 


h 
2 
Ed E,d re) E 
Ny = if o,dz = ee = yee Pe ac ae. 
nok l-p Ox Oy re) Ox 
2 
(2.5) 
h 
= Eyd,+E.d, /dv du E,d,h dv 
Ny = Oydy, = Ps agree ees.) | Ee fee ie eae 
i l~v oy Ox ) Oy 


FIGURE 2.2, Simplified wing section. 
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h 
2 Body t Epa: {0a - 
Ny =Nyx = is ds i= a re (27) 
4 ley dy Ox 


(2.8) 
aul 
se 2 
+ — E,,a0,,dz + (4, +d, + =) 
6 -v 
ci 
2 
ON, ON, 
and X= », Y=— (2.9) 
Ox oy 


where the subscripts U, L, and W refer to the upper plate, lower plate, 


and web respectively. Also, u and v are components of deflection in the 
u and v directions respectively. 


2.2 Conditions for Similarity in Mean Deflection 


Since X = ON, /dx and Y = ON, /Oy we can deduce from EQUATIONS 
2.2 the following similarity conditions: 


Ny =Ny with Ny =M, (2.10) 


EQUATION 2.4 gives 


Ne Nay Nw Ny Ny NN 


N;N 
eds, : Sa ae ae oC emeeress Qi) 
Ns N, Ns Ns v 
3 
Nay N Ne NayN 
Neg ee ee ae oe or Nav=N; (2:12) 
4n7 2 ; 
EQUATION 2.7 gives 
NgNg_,N. NENa,,N 
Wop oct Na = Ny (2.13) 
NAL NUANL 


E QUATION 2.8 gives 
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NX Ne NaN, Nav Nz NaN, Naw NW, Nau Ns 
g@2—~ tee Es Oot 2 Naw = 
N, . N,N, N, 
(2.14) 


The similarity conditions given by EquaTions 2.10, 2.11, and 2.14 can 
be summarized 


Ne Naz N,zN Np Naz 1 
“aoa or NN St (2.15) 
N,, NS v 


2.3 Local Skin Buckling of Wing Structures 


In the consideration of local buckling, we assume that there are no 
body forces or tangential forces acting in the x and y directions at the 
faces of the plate; that is, X and Y are zero in EQUATION 2.1. We then 
have, for the equilibrium equation, 


ON, ONxy F 
ox Oy i‘ 
(2.16) 
ONxy F ON, o 
ox oy 
EQUATION 2.4 simplifies to 
Eh? O*w O*w O*w 
Tos oe nee ee reas me 
ASC ey") Ox Ox? dy? oy* 
(2.17) 
‘ 0*w fe 02w oh 0?w 
p+N, + + beens: 
ox? ” dy? ”  @xdy 
where now 
du Eydy ou ov h [{ 0?w 07w 
Ne= [ody = 22 |( — +y — J -—(— py 
; l-p ox dy 2 \ dx? dy” 
(2.18) 
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or for the lower plate, we use d,, instead of dy. 


2.4 Similarity Conditions for Local Buckling 


EQUATION 2.17 gives 


NeNegN;, . NyNe |  .N; Nag, 


Vv Py $ Vv 
Ne NaN N,2 
Ed w u 
> =P oor Ny = — (2.22) 
N,4N N, 


From EQUATION 2.18 we have 


Nea, Nu NeNa NN 
Ny = ——j——_ = —— 2 ot MA =N, 2.23) 
N2Ny, N?N, 


3.0 SUMMARY AND APPLICATION OF SIMILARITY PARAMETERS 
3.1 Summary of Similarity Parameters 


It follows that to design a model is to simulate the structural behavior 
of a body in the presence of elevated temperatures except for scale fac- 
tors which are approximately one. As illustrated by similarity EQUATIONS 
1.7 and 1.8, true simulation of the deformation of thermal origin 
We, = N,N.) and that of mechanical origin (Ne, = N,/N») and can be 
only achieved with different scale factors. Of course, when one of the 
deformation effects (thermal or mechanical) is negligible compared with 
the other, the similarity in the remaining large deformation effect can be 
obtained directly through its similarity laws. 

For completeness, TABLE 3.1 summarizes the similarity laws for sim- 
ulation of aerodynamic (mechanical), elastic, inertial, and heat transfer 
effects. The complete derivation of the laws shown therein is given in 
Scipio and Teng, 1954 (Part I). 


3.2 Application of Similarity Parameters 


FicuREsS 3.1 and 3.2 show the transient aerodynamic heating phase 
‘and those regions of negligibility of either the thermal or mechanical 
effects. FicGuRES 3.1 and 3.2 are applicable to solid and wing-type 
structures, respectively. 
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TABLE 3.1 


SIMILARITY PARAMETERS FOR WIND TUNNEL TESTING OF AIRCRAFT STRUCTURES 


Phenomena investigated 


Item 


Aerodynamic Thermal Elastic Aeroelastic 
effects effects effects effects 
(1) Coefficients of (1) Surface coeffi- (1) Deflections (verti- | (1) Vibration Modes 
forces and moments cients of heat trans- cal and angular) is 
(Cr, Cy, ete.) fer and temperature (2) Critical flutter speed 
Information recovery factors, etc. (2) Strain distribu- : 
obtainable (2) Pressure distri- tions (3) Amplitudes of 
bution around the (2) Temperature dis- vibration 
model tributions (surface 2 
and internal) (4) General interaction 
Same Mach number Ng NgNaw 1 Nm =N Ne 
Same Reynolds Ne N,Ng a 2 
number ¢ Neg Ni Ny Nr = NmNz 
Approx.same Prandtl ® 2 @xiy 
number NpNLNy Ne, =NaNy/N¢ 
2 N, = 
Nd, ds N 4% 
M, = * 4 Nog =(Ner/NmNz | 
General N_P, =N ” N, fe a3 . 
similarity ipa ks NG Xs ie ier et Ny =NyNy 
criteria L page 
Ny N5=Nq eee 
N= J Na =Nz 
Nx, Nu =Ny=N5 
Same Mach number If the mechanical Nm =z 
Same Reynolds and thermal effects 3 
number are of the same order Nin =Nr 
Approx. same Prandtl : of magnitude and in Ne, =N;2 
Special number with addition, the model EI L % 
similarity and prototype materials = NO N,* 
criteria* Ng =1 are identical with Nop =(Ne1/Nm Nel 


unscaled temperatures, Ny =Nzy 
similarity cannot be 
achieved.** Np =Nzy 


* Model and prototype are constructed of identical materials, and temperature unscaled 
**See FIGURE 3.1, 
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